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Abstract 

This article surveys the mathematics of the cut and project method 
as applied to point sets, called here model sets. It covers the geometric, 
arithmetic, and analytical sides of this theory as well as diffraction and 
the connection with dynamical systems. 



Dedicated to the memory of Richard (Dick) Slansky 

The spirit of the universe is subtle and informs all life. Things live and 
die and change their forms, without knowing the root from which they come. 
Abundantly it multiplies; eternally it stands by itself. The greatest reaches 
of space do not leave its confines, and the smallest down of a bird in autumn 
awaits its power to assume form. — Chuang Tzu (tr. Lin Yutang) 



1 Introduction 

Even when reduced to its simplest form, namely that of point sets in euclidean 
space, the phenomenon of genuine quasi-periodicity appears extraordinary. 
Although it seems unfruitful to try and define the concept precisely, the 
following properties may be considered as representative: 

• discreteness 

• extensiveness 
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• finiteness of local complexity 

• repetitivity 

• diffractivity 

• aperiodicity 

• existence of exotic symmetry (optional). 

The purpose of this paper is to give an overview of the mathematics of the 
cut and project method which not only provides a very rich harvest of point 
sets (called model sets) satisfying these properties, but also provides a very 
natural way to link these ideas with many other structures in mathematics. 

A subset A of M. d is called a Delone (Delaunay) set if it is uniformly 
discrete and relatively dense. This means that there are radii r, R > so 
that each ball of radius r (resp. R) contains at most (resp. at least) one 
point of A. Although this is a fairly strong version of the first two items 
on our list, it is the most commonly used one and coincides well with the 
primitive atomic picture of a (ideally infinite) piece of material. 

The set A has finite local complexity if for each r > there are, up to 
translation, only finitely many point sets (called patches of radius r) of the 
form A PI B r (v). Here B r (v) is the ball of radius r about the point v G M. d . 
So, on each scale, there are only finitely many different patterns of points. 
This condition can be expressed topologically: A has finite local complexity 
iff the closure of A — A is discrete. It is conceivable to replace "translation" 
by "isometry" in this definition, but the theory would change considerably 



and, with the notable exception of the pinwheel tiling [|3T], little has been 
said so far on this more general situation. 

Repetitivity means loosely that any finite patch that appears, appears 
infinitely often. More precisely, given any patch of radius r there is an R so 
that within each ball of radius R, no matter its position in M. d , there is at least 
one translate of this patch. A stronger form of this requires in addition that 
each type of patch of radius r should appear with a well-defined frequency. 
The sets that we deal with here normally have this additional property (see 
Section 3). 
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From the very beginning, diffract ivity has been the hallmark of aperiodic 
order. Physically it is the most visible of its manifestations. Mathematically 
it is one of the most subtle and least visible! Very roughly we are asking 
(mathematically) that the Fourier transform of the autocorrelation density 
that arises by placing a delta peak on each point of A, should contain a part 
that looks discrete and point-like. Later on we will make this very vague 
prescription precise. The amazing thing is that in the context of model sets 
we obtain perfect diffract iveness, in the sense that the diffraction is purely 
point-like, under the fairly mild hypothesis of regularity. One of the goals of 
this survey is to show how this comes about. 

Lack of periodicity speaks for itself. Lattices and unions of cosets of 
lattices are the basis of the most prevalent forms of long-range order (crys- 
tallography). Point sets based on them satisfy all the previous properites. 
But of course, that is the trivial part of the theory! The objective is to move 
into new territory. 

Exotic symmetry usually means non-crystallographic symmetry. Although 
not mathematically essential, certainly the existence of physical structures 
with "forbidden" icosahedral symmetry was instrumental in the rapid devel- 
opment of this field. 

Even with the rather strong interpretations on the various properties 
listed above, we still do not know how to characterize sets that satisfy them. 
For an extensive discussion of these problems see |24| . However there is one 
very general method of construction which relies on controlled projection 
from a discrete group located in some auxilliary "embedding" space. In its 
original form this so-called cut and project method is based on projection 
from lattices in higher dimensional spaces. Many people have written about 
this starting, in physics, with the work of P. Kramer |21] and including the 



very useful article of Y. Meyer given in the previous edition of this School |26 



Meyer had already thought about sets formed by projection from the view 
point of harmonic analysis long before the discovery of quasi-crystals E5 



Even though it is convenient to somewhat rearrange the main components of 
his original construction, nonetheless he created a formalism which is ideal 
for creation of points sets with the desired properties of long-range aperiodic 
order. These are the model sets. 
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Some people object to the terminolgy "model set" prefering "cut and 
project set" which sounds more serious and professional. However, we prefer 
to interpret "model" as meaning exemplary and think that in terms of both 
of its priority and its greater generality the term deserves to be adopted. 

The main purpose of this article is first to give some idea of the scope 
of the relevant examples that arise as model sets (this scope surely not yet 
fully realized) and then to show how the model sets are poised between a 
number of quite different areas of mathematics. It is the satisfying way in 
which they connect many diverse parts of mathematics that makes model sets 
so intriguing and offers to the imagination so many tantalizing prospects 
for future work. For the reader interested in more on the tiling side of 
quasiperiodicity we recommend the survey paper pL which also provides 
a complementary source of some of the material presented here. 

2 Model sets 

Let us launch ourselves directly into the notion of a model set. By definition, 
a cut and project scheme consists of a collection of spaces and mappings 

K d A- W 1 x G G 

U (1) 
L 

where M. d is a real euclidean space and G is some locally compact abelian 
group, 7r 1 and tt 2 are the projection maps onto them, and L C M d x G is a 
lattice, i.e., a discrete subgroup such that the quotient group {R d x G)/L is 
compact. We assume that 71^]- is injective and that ir 2 (L) is dense in G. We 
call W 1 (resp. G) the physical (resp. internal) space. The product M. d x G is 
the embedding space. We write L = vr 1 (L). It is very convenient to define 
the mapping 

*:L-^G :^n 2 (n 1 \l 1 (x)). (2) 
Given any subset W C G, we define a corresponding set A(W) C M d by 

A(W) = {tti(x) \ xeL, tt 2 (x) G W} = {u e L \ u e W} . (3) 
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We call such a set A (or more generally any translate of such a set) a model 
set (or cut and project set) if the following condition [Wl] is fulfilled: 

Wl: W is nonempty and W = mt(W) is compact. 

For some of the deeper results we need more precise assumptions, of which 
the following are the most relevant: 

W2: The model set A is generic if the boundary dW of its window W satisfies 
is dWf]7T 2 (L) = 0. 

W3: The model set A is regular [] if dW is of (Haar) measure 0. 

The definition formalizes the notion of a point set in WL d constructed by 
projecting selected points from a lattice in some "super-space" . The points 
selected for projection are those which fall into some bounded region when 
they are projected into the complementary internal space G. The notion 
of lattice, familiar in a real space as the Z-span of a basis of that space, 
is replaced here by the more general definition that can be applied to any 
topological group. In condition [Wl] the equality could be replaced by "c" , 
but it is convenient to have this additional hypothesis since then A* = W. 

We asssume that the reader is familiar with the the most common, and 
only easily visualized examples of this, that are based on a pair of orthogonal 
axes, at irrational slopes (i.e. axes, through the origin of the standard lattice 
Z 2 in M 2 which are taken to be the physical and internal spaces, and a window 
which is an interval on the internal space e.g. [[37], [l|]. 

There are four different view points to the diagram above, which we can 
picture as follows. In the first place we have the physical space M. d and the 
point set A in it whose geometric properties are those that we wish to un- 
derstand and describe. Lattices are discrete groups inside larger continuous 
groups, and so may be thought of as arithmetic in origin. We will see later 
how in many interesting cases the arithmetic aspect is quite central. Why 
the internal side should be thought of as having to do with analysis will also 

lr rhe terminology here is not standardized. Sometimes what we call generic is called 
regular. Nonetheless, the generic situation is justifiably "generic" as we point out below. 
Our use of regular is close to the one used in |3^, |3(i| . 
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emerge in later, but an initial way to think of it is that on the internal side, 
the set of points of A appear in a totally different arrangement so that their 
closure is a very nice region of space. Finally, by definition, T := (R d x G)/L 
is a compact abelian group. In the usual situation of a real internal space, T 
is a torus, whence the notation. In any case, T has a totally natural action of 
M. d on it and it is this action that gives rise to a dynamical system. In the end 
we will have a second, related, dynamical system which plays an important 
role in questions around diffraction. 



dynamical systems side 
T 

T 

R d — R d x G — ► G (4) 

physical side j analytical side 

L 

arithmetic side 

In fact this picture can be dualized, thereby producing yet another four 
pictures! This dualization plays quite an important role in Meyer's theory 
which we touch on only most briefly here. However we will use one part of 
the dual picture. In the dual picture it is T that is the lattice and we have 

T (5) 
f 

Here we have identified the dual of the direct product M d x G with the direct 
product of the duals, and have chosen to single out the canonical projections 
as the important maps, designating them by 7Ti and TT2 respectively. 



3 Geometric side 



The geometric properties of model sets A = A(W) have been described in 
detail elsewhere (for instance [p7, 35]). We will restrict ourselves to pointing 
out a few of the most important features here. In the first place, model sets 
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are Delone sets and have the property of finite local complexity. In fact they 
satisfy a very strong form of finite local complexity: 



A — A is uniformly discrete. (6) 

A Delone set satisfying (|6|) is called a Meyer set. There are a remarkable 
number of ways describing Meyer sets |25[] ) which link them strongly 

with harmonic analysis. Though the Meyer property is considerably weaker 
than that of a model set, we nonetheless have 



Theorem 1 (W^j) Any Meyer set is a Delone subset of some model set. 

The situation regarding repetitivity is complicated by the boundary of 
the window W . A is repetitive if it is generic. If we are allowed to modify a 
model set by moving its window around then it is straightforward using the 
fact that G is a Baire space to see that the window can be moved to make the 
resulting model set generic. A proof of this can be found in J7J. Furthermore, 
in the regular case, the frequency of repetition of each patch is well-defined 
in the sense that for each finite patch P the number of occurences of the 
patch P (up to translation) per unit of volume in the ball B r (0) of radius 
r approaches a positive limit as r — ► oo. This is actually not hard to prove 
once one has established uniformity of projection (Theorem |2|). 

Lack of periodicity is automatic for model sets as long as the mapping * 
is injective. Otherwise, the kernel of * is the translation group of A. 

The comprehensive paper of Lagarias is the most extensive study to 



date of the geometry of point sets in the context of quasi-periodic structures. 



4 Arithmetic side 

Although the requirement in the definition of a model set of the existence of a 
lattice is not in itself particularly arithmetic, nonetheless the interesting and 
important examples all have strong arithmetic aspects. In the usual cases 
where the internal space is a real space, the arithmetic arises through the 
standard inner product on the embedding space and the nature of the two 
projections. 
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We will illustrate here the typical arithmetic input into the theory with 
two very different examples. 



4.1 The icosian model sets 

It was M. Baake et al |Q| who first pointed out that the root and weight 
lattices of types A 4 and D 6 could be used as the lattices for projection in 
cut and project schemes for dihedral P5 and icosahedral symmetries in 2 
and 3 dimensions. The fact that these two groups are Coxeter groups (finite 
reflection groups) and form the first two of the the series H 2 ,H 3 , if 4 of non- 
crystallographic finite Coxeter groups f\ (of which H s and are the only 
examples of rank larger than 2) suggests that H4 should also appear in this 



context. This was first pointed out in and elaborated in more detail in 
|28| . We do nothing more than outline this here. It is not necessary to know 
anything about root systems to follow this example. 

The elements of norm 1 of the usual quaternion ring H = M+Mi+Mj+RA; 
form a group isomorphic to SU(2). Since this group is a 2-fold cover of 
the orthogonal group SO (3), in particular it contains 2- fold covers of the 
icosahedral group. One such example is the following list / of 120 vectors: 

|(±1, ±1, ±1, ±1), (±1, 0, 0, 0) and all permutations, . . 

~(0, ±1, ±t', ±t) and all even permutations 

where r = (1 + y/E)/2 is the Golden ratio and ' indicates the conjugation 
map y/E i-> -y/E. 

The subring I generated by this group is called the icosian ring. Of course 
it depends on our particular choice of I, though it is straightforward to see 
that I is unique up to inner automorphisms of H. The form of the points of 
/ makes it clear that I is a Z[r]-module. We let * denote the mapping on 
I that conjugates each of the coordinates with respect to the unique Galois 
non-trivial automorphism on Z[r] (defined by sending y/E 1— > — y/E). Note 
that I* ^ I. 

2 The group H2 is the dihedral group of order 10. Usually it is fitted into the series 
12(h) of dihedral Coxeter groups, but it is also completely natural to think of it in the 
icosahedral series, as we do here. In fact, we could go a step further and include H± which 
is simply the reflection group of order 2. 
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The ring I is of rank 4 over Z[r] and rank 8 over Z. We make an explicit 
embedding of I lattice I in R 8 by the mapping x i— > (x,x*). 

This already provides the framework of a cut and project scheme: 



x 



U (8) 

i 

with the projections being given by the first and second components of (x, x*). 

Remarkably the lattice I has an entirely natural interpretation as the 
root lattice of type Eg (see for instance fill which underscores its arithmetic 



nature. This is explained in [13|, |28], |10| 



Now we wish to show that this cut and project scheme respects the sym- 
metry that is inherent in its construction. Geometrically the points of / form 
the vertices of a regular polytope P in 4-space and also form the vectors of 
a root system A 4 of type if 4. The Coxeter group is none other than the 
group of automorphisms of P (and also of A 4 ), and is in fact very easily 
described: it is the set of all (14400) maps 

x 1— > uxv ; x I— > uxv (9) 

where u,v G I. The subgroup of these transformations in which v = u^ 1 is 
obviously a copy of the icosahedral group and this subgroup stabilizes the 
3-dimensional space Ri + Rj + Rk of pure quaternions. 

These maps provide automorphisms of the rings I and, via conjugation, 
on I* too, and thus give rise to an action of I as automorphisms on the entire 
cut and project scheme. If the window W is chosen to be invariant under / 
then the resulting model set is also /-invariant. 

Restricting everything to the pure quaternions we get a new cut and 
project scheme based on the 6-dimensional root lattice D e and an icosahedral 
symmetry. Restricting further to the planes orthogonal to the 5-fold axes 
brings us back to and the related dihedral D 5 symmetry. A step further, 
and we arrive at the Fibonacci chain in 1 dimension. Thus all three families as 
well as the fundamental Fibonacci model sets fit together in this quaternionic 
model. Not only is this very pretty, it also essentially encompasses the generic 
situation for icosahedral symmetry in model sets: the only other relevant 
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lattices in 6-space are the D e weight lattice and the lattices lying between 
the root and weight lattices. For more on this see B3, ITD . 



4.2 j9-adic model sets 

Until recently, little thought had been given to the situation in which the 
internal group is something different than another real space, or at worst 
a real space crossed with a torus. However, there is a whole series of very 
natural locally compact abelian groups that are not euclidean in nature, 
namely the p-adic groups. Since these may not be familiar in this context let 
us recall the basic ideas. 

Let p be a prime number in the integers Z. Using p we can define a 
metric on the rational numbers Q, and by restriction on Z, in the following 
way. For each a G Z, we define its p- value, f p (a), as the largest exponent 
k for which p k divides a (with ^ p (0) := oo). This function is extended to 
the p-adic valuation v p : Q — > Z by u p (a/b) := v p (a) — v p (b) for all rational 
numbers a/b. We now define the "distance" between two rational numbers 
x, y as d(x, y) = p~ v vky- x ) . 

It is not hard to see that this does define a metric on Q, in which closeness 
to is equivalent to high divisibility by the prime p. The completion of 
the rationals under this topology is the field of p-adic numbers Q p and the 
completion of the subring Z is the subring of p-adic integers, Z p . Each p- 
adic integer can be given the more concrete representation as a series in 
the form Yl^=o a nP n where the a n are integers in the range < a n < p. 
Note that convergence here is automatic, even though there are infinitely 
many terms in the sum, because of the nature of the p-adic topology. The 
topologies defined by such metrics have other counter-intuitive properties. 
For example, for each non-negative integer k, the set p k ■ Z p , of elements of 
Z p divisible by p k , is the ball of radius p~ k and is clopen, i.e. both open and 
closed, as too are all its cosets, a + p k ■ Z p . 

Seen as a topological space, Z p is both compact and totally disconnected 
(but not discrete). In particular, Q p and Z p are locally compact abelian 
groups under addition. Thus, we can use Z p to construct interesting cut and 
project schemes for M. d simply by taking G := (Z p ) d and L = Z d embedded 
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diagonally into M. d x Z p d (based on the natural embedding of Z in Z p ). For 
more on p-adic numbers and other totally disconnected groups, the reader 
may consult [215, [|. 

In it was shown that a number of interesting substitution systems and 
tilings can be interpreted in a p-adic setting, including the well-known chair 
tiling. Rather than repeat these, let us give a different example, mentioned 
in |7| but not elaborated upon. One of the earliest classes of aperiodic tilings 
to be discovered was the class of Raphael Robinson's square tilings |33| . The 
title of his paper recalls that the mathematical interest in aperiodic structures 
had a totally different (and earlier) origin than the physical one, namely the 
interest in decidability problems in the tiling of the plane with tiles of finitely 
many different types. The Robinson tilings are tilings of the plane by equally 
sized squares, in the usual fashion, with the twist that the square tiles come 
in 6 types (up to rotational and reflectional symmetries), distinguished by 
the markings of their edges, and the tiling is required to respect these edge 
markings by having the edges of adjacent tiles properly matched. Pictures of 
the tiles may be found in [33] and [TjJ. What is important for our discussion 
is that there is another set of markings by lines of these tiles and the correct 
tilings are those for which these lines arrange themselves into a pattern of 
squares of increasing scales 1,2,4,8, . . . (see Fig.l). This picture is the one 
that Robinson used to prove the aperiodicity, for evidently no translation 
can map the squares of all scales onto themselves simultaneously. The same 
idea was used by Penrose in his recent hexagonal tiling pOi . 



Now the point is that the centres of the tiles of each of the six types 
form a model set based on an internal space which is 2-adic. Very briefly the 
argument is as follows. 

Starting with Fig. 1 and ignoring the actual square tiles themselves we 
have patterns of interlocking squares of increasing scale. Let's call these 
the pattern- squares to distinguish them from the actual tiling squares. The 
vertices of the smallest scale pattern-squares (order 1) form the vertices of a 
lattice V once one of them, say c\ = (0, 0), has been chosen as the origin. It is 
convenient to introduce the larger lattice L — which we can identify with 
Z 2 . The locations of the vertices of the increasingly scaled pattern-squares 
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Figure V. Pattern of increasing squares in a Robinson square tiling 



determine two sequences, {a k } and {j3 k }, composed out of the two numbers 
{±1} as follows: the vertices of the pattern-squares of order k (side-length 
2 k ) are the points of a coset 

L k :=c k + 2 k L. (10) 

where 

c k = (a + ai2 + . . . , +a fc - 2 2 fc ~ 2 , ft + ft2 + . . . , +/3 fc _ 2 2 fe - 2 ) (11) 

for k = 2, 3, . . .. Conversely, given two sequences a, /3 of ±l's we can use (|TTJ) 
and ([ToD to define the vertices of a suitable pattern of squares. 

The pattern-squares themselves are determined by the condition that the 
points of L k are the centres of the squares of order k — 1 for all k > 1. 
This then establishes a coordinatization of the pattern-squares. Both the 
actual pattern of the pattern-squares and their coordinatization depend on 
the choice of our two sequences (more below). 

Now looking again at the tiling squares, we distinguish the 6 types of 
tiles according to their location in the pattern-squares. We list these here 
together with a description of the coordinates of the centres of their tiles: 

(1) the "corner tiles" of the squares of order 1; 
coordinates L x = 2L; density i; 
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(2) the "corner tiles" of all squares of all higher orders; 
coordinates Ufe>2 c k + 2 fc £; density ^ 

(3) "cross tiles" where edges of two different orders of pattern squares meet 
(actually these orders always differ by exactly 1); 

coordinates: 

(U> 3 Ok + (±2^ 2 , ±2 fc - 3 ) + 2 k L) U (U> 3 c k + {±2 k -\ ±2 fe - 2 ) + 2 k L) ; 
density |; 

(4) "edge " squares, which contain part of a single edge of a pattern-square, 
except those in which are exactly in the middle of an edge; density |; 

(5) "edge " squares, which contain part of a single edge of a pattern-square 
and which are exactly in the middle of an edge; density |; 

(6) blank tiles, with no part of any edge in them; density |. 

Observe that each of these sets is a countable union Wj of cosets of the 
form a + 2 h L. Let us replace each of these by the corresponding 2-adic clopen 
set a + 2 h -Z2- In this way we get 6 open sets Wj, j — 1, . . . , 6, whose closures, 
being closed subsets of the compact group Z 2 , are compact with non-empty 
interiors. Finally we can describe the centres of the squares of type j as 

{x G Z 2 | x G Wj} (12) 

which, unlikely as it appears, is a model set under the scheme 

R 2 R 2 xZ 2 Q 

U (13) 
Z 2 

where Z 2 is embedded into M 2 x Z 2 . diagonally: x i— > (x,x). 

The entire tiling is determined by the vertices of the various squares and 
hence by the window 

oo 

W:=\Jc k + 2 fc Z 2 . (14) 
k=i 

Evidently W is an open subset of Z 2 .. Let w G Then for each m G Z + , 

w + 2 m Z 2 . meets W, so w = c k mod 2 min ^ m ' fe ^ Z% for some k = k(m). If k < m 
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then w G c k + 2 k l? 2 C W. Thus k > m and w = c fc ( m )mod 2 m Z 2 ,. It follows 
that w is the limit of some subsequence of the {c^}. Since the entire sequence 
evidently converges (in the p-adic topology, of course!) to some c = (a, b), 
where a = a fc2 fc and similarly for b, we see that w = c and so dW = {c}. 
Thus the model set of all vertices is regular, and even generic provided that 

ci z 2 . 

More generally, one may expect this p-adic topologies to arise whenever 
there is a self-similarity 9 : L — > L for which 6(L) C L, but 0(L) ^ L. 

In [0 we also see the appearance of mixed p-adic and real spaces as the 
internal spaces. Beyond these types we are not aware of any interesting 
examples, though they may well exist. 



5 Analytic side 

The transition from an inherently discete picture on the physical side to 
something inherently far more continuous on the internal side is made via 
H. Weyl's theory of uniform distribution. 

Let us assume that we have a model set A. Now consider the following 
question. Suppose that we take a ball -Br(O) of radius R about the origin in 
M. d and look at A# := A fl Br(0). Then we can ask how A* R is distributed 
over W. We say that the sets A* R are uniformly distributed if for each open 
set U C W we have 

cardfApfW) . . TJT , . . 

where fi is Haar measure on G. 

Theorem 2 j^h] , |7^ / If A is regular then the sets A* R are uniformly distributed 
over W . 

Let /* : G — > C be any function. We define / : L — > C by f(x) = 
f*(x*). If /* is supported on the window W then evidently / is supported 
on the model set A. If /* is continuous (which is the case of interest) then 
this is iff. 
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Theorem 3 (Weyl) j^j If A is regular and f* is continuous then 

lim i— V f(x) = — TTfTTT / f*(uW(u) (16) 

R-+00 card(A fl ) v ; vol(W) 7vk 

Since has boundary of measure zero, it is not necessary to insist that 
/* (which is supported on W) be continuous on all of internal space, only on 
the window W. 

In this way discrete averaging on the model set is transformed into inte- 
gration on the window. This process was used in [|j] to determine the exis- 
tence of invariant measures on internal space in the presence of self-similarity 
on the quasi- crystal. We briefly explain this. We assume here that internal 
space is R n for some n. 

A self-similarity of A is an affine linear mapping t = tq v 

tq,v : x^Qx + v (17) 

on M. d that maps A into itself, where Q is a (linear) similarity and v G M. d . 
Thus Q = qR, i.e. it is made up of an orthogonal transformation R and an 
inflation factor q. 

Let tq v be a self-similarity of A. Since A is uniformly discrete, we must 
have \q\ > 1. We will assume \q\ > 1 and that QL = L. We are interested in 
the entire set of affine inflations with the same similarity factor Q. 

Note that Q naturally gives rise to an automorphism Q of the lattice L, 
i.e. an element of GLz(L), and a linear mapping Q* of MJ 1 that maps W into 
itself. From the arithmetic nature of Q we deduce that the eigenvalues of 
Q and Q* are algebraic integers and from the compactness of W that Q* is 
contractive. 

Define 

W Q := {u e W 1 | Q*W + u C W} , (18) 

We say that Q is compatible with A if int(Wg) ^ 0. Assuming that this is 
the case (not a strong assumption) , then the set Tq of affine inflations with 
the same similarity Q is the set of mappings tg : x \— > Qx + v, where v runs 
through the set 

T = T Q := {v G L \ v* E W Q } . (19) 
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Theorem 4 If Q is a self- similarity and the above assumptions on Q apply 
then there is a unique absolutely continuous positive measure // on internal 
space, supported on W , satisfying: 

• »(W) = 1; 

• fi is invariant in the sense that, if we define t* • /ij by t* ■ n(Y) = 
ti((fi)- l (y)), then 

a = Urn— — — V t* r u f . (20) 

7TK V " v£TDB s (0) 

The similarity of this measure to Hutchison measures in the context of iter- 
ated function systems is not coincidental. In fact, if we restrict to the ball 
B s (0) then the {£*} form a finite set of contractions which is indeed an iter- 
ated function system. For more on this and invariant density functions on 
model sets see [|]. 

The type of limit averaging involved here is a very natural one from the 
point of view of physical situations, representing the transition from the world 
of sets finite in extent to the ideal world of infinitely extended point sets. 

Although no one to our knowledge has made any use of it, it is interesting 
to use Weyl's theorem to transfer the structure of L 2 (W) to a space of similar 
objects on A. Namely, the space of continuous functions on W leads to a 
space 

C = C(A) (21) 

of funtions on A via the mapping *. Then the usual inner product (f,g) = 
J w f*{u)g*{u)du = (f,g)w defines an inner product on C and we can com- 
plete this space in order to get a Hilbert space C isomorphic to L 2 (W). Of 
course the elements of C can no longer be interpreted as functions on A since 
functions on A that differ by a function whose absolute square has limit 
average sum equal to are identified. 

6 Dynamical systems side 

So far we have looked at one model set in isolation. Now we move on to 
consider families of model sets. We start with a number of definitions and 
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results. All of these may be found in the paper of M. Schlottmann on 
which we have relied heavily here. Many are well-known in the context of 
tilings for which a recent reference with a good bibliography is |33|. In this 



section all point sets under discussion are assumed to be Delone sets in M d . 

Two Delone sets S, S' in M. d are locally isomorphic (or some people say 
locally indistinguishable) if, up to translations, every patch of either of them 
occurs in the other. Thus on any finite scale, up to translation, the two sets 
are indistinguishable. Given a Delone set S C M. d we can look at its local 
isomorphism class (LI class) LI(S'), namely all point sets locally isomorphic 
to it. 

We denote by X(r) the set of all Delone sets of M d for which the minimum 
separation between distinct points is at least r. We assume in the rest of this 
section that r > has been fixed. 

We define a Hausdorff topology on X{r) as follows: two sets S, S' £ X{r) 
are "close" if for some large compact set K C M. d and some small e we have 

(v + S) n K = S' n K (22) 

for some v £ lR r with \v\ < e. More precisely we define a uniformity U on 
X(r) using as the sets U (K, e) of uniformity the set of pairs (S, S') satisfying 

m 



Theorem 5 [3%., W\, W6 J With respect to this topology X{r) is a complete 
Hausdorff space. 

Let S £ X{r). Then IR d acts on LI(S') by translation and in particular 
the entire orbit [S] of S lies in LI(S'). This action is continuous and hence 
extends also to an action on the closure Ll(S) of LI(5). The relationship 
between orbits and LI classes can be summed up by 



S £ [S] C U(S) C [S] = U(S). (23) 

The second inclusion follows easily from the definitions. The inclusions may, 
according to the situation, be strict or actual equalities. For a lattice there 
is only one orbit in its LI class. For general model sets the situation is very 
different, as we shall see. 
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Recall that a Delone set S is said to be of finite local complexity if the 
closure of S — S is discrete. Finite local complexity is a property that is 
inherited by whole LI classes. 

Theorem 6 An LI class is pre- compact (i.e its completion is compact) iff it 
has finite local complexity. 

Thus, if S is a Delone set of finite local complexity we obtain a dynamical 
system T>(S): 

R d x[S] — >\S\. (24) 

In the sequel we will use the symbols like T>(S) to denote both the dynamical 
system itself and the corresponding defining space [S]. 

Theorem 7 Let S be a Delone set of finite local complexity. The following 
are equivalent: 

(i) S is repetitive; 

(ii) [S] = LI(S') is closed; 

(Hi) The dynamical system V(S) is minimal. 

We recall that minimal means that every M d orbit is dense. 

Since generic model sets are repetitive, this leads to a very nice result: 

Theorem 8 Let A be a generic model set. Then its LI class LI(A) is a 
compact Hausdorff space and under the action of translation under M. d it 
becomes a minimal dynamical system, P(A). 

This is the first of the dynamical systems that we wish to consider. Its 
rather abstract form is better understood by relating it to a more accessible 
dynamical system. 

To this end, let A = A(W) = {x G L \ x* G W} be a model set. Each 
element (u, v) of the group ~R d x G can be used to form a new model set 

A(W, u,v) := u + {x G L | x* G -v + W}. (25) 
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If (u, v) G L then v = u* and we can rewrite this as {u+x G L | (u+x)* G W}, 
which is just A again. Thus we get a whole family of model sets parametrized 
by T := (Mr x G)/L with M. d x G acting on it. This is the second dynamical 
system. Its points correspond to the model sets A(W,u,v). This is the so- 
called torus parametrization introduced by Baake et al. in ||. We use the 
same terminology in the more general context here, although in general T is 
not a torus! 

The action of l d on T, (x, y + L) t— > x + y + L, is a faithful transcription 
of the operation of translation in physical space, so the orbits of M d on T 
correspond to model sets that differ only by translation. The action of G on 
T corresponds to translating the window around. 

Theorem 9 Let A be a generic model set. Then 

K d xT — > T (26) 

is a minimal uniquely ergodic dynamical system T> tOT . The unique invariant 
probability measure is normalized Haar measure. The set of points of T> tOT 
corresponding to generic model sets is dense and indeed the set of points 
corresponding to the non-generic model sets is of the first category. 

It is noteworthy that this dynamical system is independent of W but the 
actual parametrization of model sets is clearly dependent on it. 

So now given a generic model set A, there are two dynamical systems 
for the group M. d , one T>(A) coming from the closure of the orbit of A under 
action of M> d and another V tOT coming from the torus parametrization. Not 
surprisingly they are related, but rather surprisingly this relation is somewhat 
subtle. All the elements of T>(A) are, by definition, in the same LI class. The 
same is not the case for the model sets parametrized by T> tOT . Indeed, A is 
generic, but translating the window around is bound to produce model sets 
that are not generic. These non-generic model sets are not locally isomorphic 
to the regular ones, because they have certain special local configurations of 
points that are related to the boundaries of their windows. 



Theorem 10 |f| Let A be a generic model set. Then there is a continuous 



surjective mapping 

(3 : V(A) — > Aor (27) 
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which is M. d -equivariant and which maps A onto the point of the torus. 
Furthermore, for each of the points of T> tor which parametrize generic model 
sets, the preimage in P(A) consists of a unique point. 

This mapping comes about as follows: Let A' e T>(A). First suppose that 
A' C L. Then it is not hard to see that PLeA'C^ ~ x *) * s a single point, 
call it b(A'). Furthermore, for all u 6 L, b(A' — u) = b(A') + u*. Now for 
arbitrary A' G T>(A), we can always find v e M d with A' — v C L. This t> is 
nothing like unique but it follows from what we have just said that the pair 
P(A') := (v, b(A' — v)) is unique mod L, and this is the mapping that we 
require. 

Using these facts it can be established that 
Theorem 11 |36| Assume that A is a regular and generic model set. Then 



T>(A) is uniquely ergodic and furthermore L 2 {V(A)) and L 2 (V tor ) are isomet- 
rically isomorphic as M. d -spaces. 

The importance of this is that it shows that from the spectrum of T> tOT 
being discrete, which it surely is since T a compact abelian group, it follows 
that the spectrum of X'(A) is also discrete. It is from this that the pure point 
diffractivity of A can be deduced. In the final section we briefly describe how 
this happens. 

7 Diffraction 

The theoretical framework for the discussion of diffraction has been very well 
described in several places. The two papers of A. Hof [O, Il8|] are standards 



and there are also good descriptions in |15|, |J. Here we just quickly formulate 
the definitions. 

Let A be a regular model set and define the (tempered) distribution 

Sa:=^2S x , (28) 

where 5 X is the Dirac measure at x. For each s > we calculate the auto- 
correlation of 5a restricted to the ball of radius s: 

^AnB s (o) * ^AnB s (o) = /2 ^ x ~ y ' 

x,yeAnB s (o) 
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where, as usual, the over-tilde indicates changing the sign of the argument. 
The limit as s goes to infinity of the volume-averaged auto-correlation of this 
measure, which exists for model sets, is the auto- correlation measure of A 
(its so-called Patterson function): 

This limit, taken in the vague topology, converges to a tempered distribution 
(i.e. this limit exists when taken against rapidly decreasing test functions). 
Its Fourier transform is a positive measure 7 (a result of Bochner's theorem 
applied to the positive definite distribution 7) which is the diffraction pattern 
of A. The measure decomposes into a point part and a continuous part. The 
point part of this measure is the Bragg spectrum of A. The model set has 
pure point spectrum if the continuous part is the trivial 0-measure. 

The complexity of the definition makes it hard to discover the nature 
of the diffraction pattern, in particular whether or not we have pure-point 
diffraction or not. One approach has been to use the ergodic theory outlined 
above, and indeed it is able to give the main result: 



Theorem 12 |36| Any regular model set has pure point spectrum. Further- 
more this spectrum is supported on the projection into Fourier space on the 
physical side of the dual of the compact group T @ ; i.e. it has the form 

7 = ^2w(k)8z lik ) (31) 
kef 



The proof of this is based on an idea of Dworkin [|12|]. The argument is 



spelled out in |L8| and we repeat it here since otherwise it is difficult to see 



the connection between dynamical systems and diffraction. 

We can assume that A is generic since translation of the window does not 
alter the qualitative nature of the diffraction. The next step is to replace 
5\ by a smooth approximation to it. To this end, let b : ~R d — > IR> be a 
smooth function whose support is contained in the ball B r (0) of radius r, 
where B 2r {0) n (A - A) = {0}. Define a function ip : £>(A) — > R by 

ip{A') = [ b{-u)6 A ,{u)du, (32) 
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which is continuous on T>(A). The action of M d on [A] gives rise to a corre- 
sponding action x i— > T x on the space of functions on the orbit [A] of A under 
translation. 

For each x G M d we have 

T x ijj(A) = / 6(— w)5_ ;r+ A(u)(iM = / b(— u)5^{x + u)du = b * 5\(x) , (33) 

which shows that the function a^ b > : M. d — > C defined by x i— > T 3; ('0)(A) is 
obtained by centering a copy of 6 at each point of A. 
Now consider the autocorrelation of a^: 



7 W(x) = lim^oo vol( i s(0)) f BM T x+y {^){A)T y {^){A)dy 

(34) 

The main point here is the use of the Birkhoff ergodic theorem and the ergod- 
icity of the action of M d on P(A) to replace the integral over R d by an integral 
over V(A). Note that the uniqueness of ergodicity and the continuity of ip is 
needed here to guarantee the statement for all x rather than a.e. (IB, Sec. 
3.2). 



In view of Theorem [11], we have a Fourier expansion of ip in terms of the 
eigenf unctions for the action of lR d : ip = a A0A, where cj)\ is the eigen- 
function for the character x \— > e 2mX ' x on IR d . Thus (T x tfj,tp) = ^ \a\\ 2 e 2 ' x 
and taking Fourier transforms we have 

7 (6 > = 5>a| 2 *a (35) 

which is a pure point measure on M. d . 

On the other hand we know that = b * Sa whose autocorrelation can 
be calculated directly as 6*6*7 an d so 7^ = |6| 2 7. Finally, taking a sequence 
of bump functions {b} converging in the vague topology to 5 we obtain the 
required pure-point nature of the diffraction pattern. 

Theorem |I2| is qualitative in nature. The quantitative counterpart is this: 



Theorem 13 fi2by Let k 6 T and let \ denote the characteristic (or indica- 
tor) junction of W . Then w(k) = \x(~n 2 (k)) /vol(W)\ 2 . 
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There are a number of variations on this theme that are worthwhile men- 
tioning. First we may imagine replacing the simple sum (^Sj) by a weighted 
sum 

#a := y^ y m(x) 5 X , (36) 

xeA 

where u : A — > C is some function. 

Theorem 14 If A is a regular model set and if uj G C(A) (see ft2H{)) then the 
weighted point distribution is pure point diffractive. 

Next we consider the case that our points of A are considered stochasti- 
cally: 

^stochastic := ^ V{ X ) $x 5 (37) 

where the r](x) form a collection of independent identically distributed ran- 
dom variables that take the values 1, (indicating occupancy or not of the 
respective model set sites) with the probability of occupancy being p: 

Theorem 15 Let A be a regular model set and let rj be as above with 
the mean and second moment equal to mi and m 2 respectively. Then the 
autocorrelation of A and that of its stochastic version are, with probability 
one, related by 

7stochastic = (^i) 2 7 + d (m 2 - (rai) 2 )5 (38) 
with Fourier transforms 

7stochastic = (mi) 2 7 + d (m 2 - (^i) 2 ) , (39) 
where d is the density per unit volume of the A. 

Thus the pure point nature of A is affected by at most the addition of a 
constant continuous background. More on this stochastic approach may be 
found in JJ. 

Yet a different variation is to allow the points of A to be moved in some 
regular way. 
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Theorem 16 j\T^ Let A be a regular model set and let f : x i— > (fi(x), . . . , fd{x)) 
be some mapping of A into M. d , where each f\ G C(A) (see Then the 

set Af := {x + f(x) \ x G A} zs jmre pomt diffractive. 

These variations can be combined in the obvious ways. 

The problems of determining which point sets are pure point diffractive 
is a fascinating and challenging one which is still wide open. The examples 
above show how much model sets can be modified without serious damage 
to their diffractive properties. Obviously adding or removing points whose 
average density is also does not alter the diffraction. But there are point 
sets that are even more remote that are diffractive. One example is the set of 
visible points of a lattice. Given a lattice L in IR d its visible points are those 
points x G L, x ^ 0, satisfying Qx n L — Zx. What is interesting about the 
visible points is that they do not form a Delone set (they are not relatively 
dense in R d ). In fact, for each r > 0, the set of holes of radius exceeding r 
has positive density. However, 

Theorem 17 ^ The set of visible points of any lattice of rank at least 2 is 
pure point diffractive. 



7.1 Comment 

In this paper we have considered point sets in M d that are constructed through 
the method of projection from an embedding group and a lattice. We have 
assumed that the embedding group is of the form tfxG where G is a lo- 
cally compact abelian group, which, in view of Theorem is fairly natural. 
However, it is possible to study model sets in the situation where the 'phys- 
ical space' has been replaced by an arbitrary locally compact abelian group 
without losing many of the most interesting properties. In particular the 



diffraction results of Theorem 12 have formulations in this generality fl3q 
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